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摘 Xe: 水 文 模型 参数 的 敏感 性 分 析 、 优 化 和 验证 对 提高 模型 率 定 效率 和 计算 精度 具有 重要 意义 。 为 探讨 Laio 土壤 水 
分 动态 随机 模型 (Laio 模型 ) 各 参数 在 垄沟 集 雨 系统 的 敏感 性 ， 同 时 确定 参数 优化 和 模型 验证 的 最 佳 方案 ,本 文 结合 多 因 
素 敏 感性 分 析 法 以 及 改进 单纯 形 法 (ISM)、 粒 子 群 优化 算法 (PSO) 和 混合 粒子 群 优化 算法 (HPSO), 利用 中 国 气象 局 定西 干 
时 气象 与 生态 环境 试验 基地 2012—2013 年 垄沟 集 雨 燕麦 生长 季 降 雨 、 径 流 和 土壤 水 分 等 实测 数据 ， 对 垄沟 集 雨 系 统 
Laio 模型 的 13 个 参数 进行 敏感 性 分 析 、 优 化 和 验证 。 结 果 表 明 ,平均 降水 量 w MAERA ww 对 土壤 水 分 概率 密度 函数 
p(s) 最 敏感 ， 且 p(s) 对 参数 a 的 敏感 性 在 低 土 壤 含水 率 下 更 明显 ， 而 对 参数 ww 的 敏感 性 在 高 土壤 含水 率 下 更 明显 ; 3 种 
算法 ISM、PSO 和 HPSO) 的 优化 参数 值 均 能 对 垄沟 集 雨 系统 土壤 水 分 概率 密度 函数 进行 较 好 模拟 ,其 峰值 (CPV)、 峰 值 
位 置 (PP) 和 95% 置 信 区 间 (CI95%) 实 测 值 与 模拟 值 的 相对 误差 均 小 于 10%, CM 指数 均 大 于 0.5; 同时 , HPSO 算法 优化 参 
数 的 模拟 效果 和 收敛 速度 均 显 著 优 于 PSO 算法 和 JSM 算法 ， 且 能 较 显 著 克 服 ISM 算法 和 PSO 算法 存在 的 缺陷 。 故 HPSO 
算法 可 作为 垄沟 集 雨 系统 土壤 水 分 动态 随机 模型 参数 优化 的 待 选 方案 。 该 研究 可 为 Laio 模型 的 区 域 应 用 和 模型 调 参 提 
供 科学 指导 依据 。 

关键 词 : 垄沟 集 雨 系统 ; 土壤 水 分 动态 ; Laio 土 壤 水 分 动态 随机 模型 ; 敏感 性 分 析 ; 模型 参数 优化 
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Abstract: Ecosystem dynamics in arid and semiarid climates are strongly dependent on the soil water availability which, in turn, is 
the result of a number of complex and mutually interacting hydrological processes. Furthermore, soil moisture dynamics is an 
indispensable part of the research on the interactions and feedbacks between hydrological processes and terrestrial ecosystem 
processes, which is the result of non-linear interactions among a series of hydrological, climatic, and ecological processes. 
Consequently, soil moisture content and its dynamics needs to be studied by stochastic eco-hydrological models, which can 
reasonably describe the characteristics of soil moisture dynamics including the pulse, erratic, and random processes. Moreover, 
eco-hydrological model calibration and validation are essential for model evaluation and application. It is important for model 
application to accurately estimate the values of eco-hydrological model parameters and to further improve the capacity of model 
prediction. In this study, to establish an effective method system for sensitivity analysis, calibration and validation of Laio soil 
moisture dynamic stochastic model (Laio model) parameters in a ridge-furrow rainwater harvesting system in semi-arid area, a 
field experiment with a randomized complete block design was conducted during the 2012 and 2013 oats growing seasons at the 
Dingxi Arid Meteorology and Ecological Environment Experimental Station. The experiment was designed to investigate the 


parameters' sensitivity and determine the optimal plan for the parameters' optimization of Laio model in different mulching 


materials (common plastic film, biodegradable mulching film and manual compacted soil) and different ridge-furrow ratios (60 : 
30, 60 : 45, and 60 : 60 [cm : cm]) of ridge-furrow rainwater harvesting system based on the multi-factor sensitivity analysis 
method and the improved simplex method (ISM), the particle swarm optimization algorithm (PSO) and the hybrid particle swarm 
optimization algorithm (HPSO). The results indicated that: (1) the mean amount of precipitation per rainfall event (a) and the 
degree of soil saturation at wilting point (sw) were the most sensitive parameters for the probabilistic density function of soil 
moisture (p(s)) in different experimental treatments, and the sensitivity of p(s) to parameter a and swis more obvious under soil 
moisture content and high soil moisture content, respectively; (2) A good agreement between the result of modeling by the 
optimized parameters of Laio model using three optimization algorithms (ISM, PSO and HPSO) and observation was found based 
on the curve shape of p(s) which included the curve peak value (CPV), the position of the peak (PP), the confidence interval of 
95% (CI95%) and the consistency measure (CM), indicating the optimized parameters of Laio model by the ISM, PSO and HPSO 
is able to estimate the p(s) in ridge-furrow rainwater harvesting system; and (3) HPSO not only improves the global optimization 
performance, but also quickens the convergence speed and obtains robust results with good quality, which indicates HPSO is an 
effective optimization method for Laio model calibration and validation in ridge-furrow rainwater harvesting system. This study 
could improve the efficiency of model parameters calibration, upgrade the accuracy of model simulation result, and provide 
guidance for the application research of Laio model in ridge-furrow rainwater harvesting system. 

Keywords: ridge-furrow rainwater harvesting; soil moisture dynamic; Laio soil moisture dynamic stochastic model; sensitivity 


analysis; model parameters optimization 


水 资源 短缺 是 全 球 旱 地 农业 面临 的 共同 难题 。 如 何 充分 利用 天 然 降雨 资源 、 确 保 粮 食 生产 安全 和 维 
持 农田 生态 系统 稳定 是 实现 旱 农 区 “农业 -生态 -经 济 ” 耦 合 系统 协调 、 持 续 发 展 的 关键 册 。 同 时 ， 研 发 和 推 
广 高 效 、 低 廉 、 环 保 的 旱 作 栽培 技术 对 提高 旱 农 区 作物 产量 及 水 分 利用 效率 具有 重要 的 意义 趾 。 垄 沟 集 
雨 系统 利用 田间 起 垄 、 沟 玖 相间 、 垄 面 产 流 、 沟 内 高 效 集 雨 ， 并 依靠 增 温 、 抑 燕 \ 保 土 等 生理 生态 效应 , 已 
成 为 水 分 缺乏 的 半 干 旱 区 农田 生态 系统 一 项 重要 的 集 水 节 汉 措施 BHA。 其 在 绥 解 旱 农 区 人 口 急剧 增长 、 粮 
食 日 益 紧缺 与 农业 生态 恶化 的 矛盾 中 发 挥 着 至 关 重 要 的 作用 。 

土壤 水 分 是 半 干 旱 农田 生态 系统 植物 水 分 的 主要 来 源 , 是 养分 循环 和 流动 的 载体 , 在 土壤 -植被 -大 

气 系统 物质 和 能 量 转化 中 起 着 核心 和 纽带 的 重要 作用 名。 认识 半 干 旱 农田 生态 系统 与 土壤 水 分 关系 和 相 
互 作用 机 理 ， 对 揭示 农田 生态 系统 稳定 性 及 其 水 土 关 键 要 素 的 变化 过 程 具有 重要 意义 四。 由 于 影响 土壤 水 
分 动态 的 各 因素 (降水 、 蒸 散发 、 土 壤 异 质 性 、 地 形 等 ) 具 有 随机 性 ,特别 是 降雨 事件 发 生 及 降雨 量 分 布 的 
随机 性 ,决定 了 土壤 水 分 动态 模型 只 有 以 概率 形式 描述 才 具 有 实际 意义 中 。 自 Eagleson 等 四 首次 将 随机 理 
念 纳入 土壤 水 量 平衡 方程 ,， 随后 Mily, Rodriguez-Iturbe 等 00、Laio 等 0、Porporato $5U2!, Pan 等 (13 
先后 对 不 同时 空 尺度 的 土壤 水 分 动态 建立 了 随机 数学 模型 ， 并 得 到 广泛 的 应 用 。Laio 土壤 水 分 动态 随机 
模型 (Laio 模型 ) 在 蒸 散 发 项 上 引进 了 2 个 土壤 水 分 临界 值 ( 凋 蓉 系 数 和 吸湿 系数 )， 其 能 对 干旱 半 干 旱 区 农 
田 生 态 系统 土壤 水 分 动态 进行 更 真实 地 描述 , 可 为 旱 农 区 土壤 水 分 的 有 效 利 用 与 管理 提供 理论 指导 。 
由 于 半 干 旱 区 垄沟 集 雨 系统 土壤 水 分 与 多 种 物理 、 化 学 和 生物 过 程 以 及 降雨 、 径 流 、 菩 散发、 土壤 
特性 、 微 地 形 及 履 盖 材料 等 密切 相关 , 长 期 处 于 动态 变化 状态 ， 且 变化 过 程 比 较 复 杂 04469， 所 以 利用 Laio 
模型 对 该 系统 土壤 水 分 动态 进行 模拟 和 研究 是 很 有 必要 的 。Laio 模型 共 涉 及 土壤 、 植 被 和 气候 等 13 个 参 
数 ， 部 分 参数 很 难 通过 观测 直接 获取 ,， 且 参数 值 存在 极 大 不 确定 性 475。 故 在 模型 应 用 前 需 考 虑 模型 “本 
地 化 * 和 “区 域 化 ”问题 , 即 需要 对 模型 参数 进行 敏感 性 分 析 和 优化 。 目 前 针对 Laio 模型 参数 的 敏感 性 及 获 
取 方 法 已 有 一 些 研究 ， 如 姚 淑 起 等 中 在 科尔沁 沙 地 对 Laio 模型 参数 的 敏感 性 进行 了 分 析 ， 并 将 参数 按 敏 
感性 强 弱 分 为 了 3 类 ; Miller 等 外 对 Laio 模型 参数 的 获取 及 其 敏感 性 分 析 发 现 ， 最 大 蒸 散 量 Ema 和 水 分 胁 
迫 点 8 最 难 获取 , 且 其 敏感 性 也 最 高 ; 任 庆 福 等 0 利用 PSO 算法 对 太行 山 山 前 平原 典型 井 灌 农 区 Laio 模 
型 参数 进行 了 优化 , 发 现 经 率 定 的 参数 能 更 好 地 模拟 作物 生长 期 土壤 水 分 的 随机 变化 特征 。 然 而 , 在 半 干 
早 区 对 不 同 含水 率 条 件 下 垄沟 集 雨 系统 Laio 模型 参数 的 敏感 性 ， 以 及 各 种 优化 算法 对 模型 参数 优化 的 适 
应 度 和 有 效 性 尚 缺 乏 系 统 研究 。 鉴 此 , 建立 方便 可 行 的 模型 参数 优化 、 敏 感性 分 析 和 有 效 性 验证 的 方法 
体系 对 提高 参数 率 定 效率 , 控制 模型 计算 误差 和 拓宽 模型 应 用 领域 具有 极其 重要 的 意义 。 

本 文 利 用 中 国 气 象 局 定西 干旱 气象 与 生态 环境 试验 基地 2012—2013 年 垄沟 集 雨 燕麦 生长 季 降 雨 、 径 
流 和 土壤 水 分 等 实测 数据 资料 , 采用 多 因素 敏感 性 分 析 法 ， 对 半 干 旱 区 垄沟 集 雨 系统 Laio 土壤 水 分 动态 
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随机 模型 参数 的 敏感 性 i 


期 建立 垄沟 集 雨 系统 Laio 模型 参数 敏感 性 分 析 、 优化 和 验 订 


域 应 用 提供 科学 理论 依据 。 
1 材料 与 方法 
1.1 ”田间 试验 和 数据 测定 


行 分 析 和 分 类 。 同 时 ,基于 ISM, PSO 和 HPSO $37 
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E ， 对 垄沟 集 雨 系统 Laio 模 
型 的 13 个 参数 进行 优化 和 优选 ,并 利用 实测 数据 资料 对 3 种 算法 优化 参数 的 有 效 性 进行 验证 和 评价 ， 以 


E 的 有 效 方 法 体系 , 为 Laio 模型 参数 校正 和 区 


试验 于 2012 一 2013 年 在 中 国 气 象 局 兰州 干旱 气象 研究 所 定 


N, 104°35' E, 海拔 1 896.7 m) 进 行 , 该 基地 属 


mm。 降 雨 在 年 内 分 布 极 不 均匀 ,7~10 月 降雨 量 


mm. 试验 地 地 势 平 坦 ,土壤 为 重 壤土 , 0~100 


盖 作 为 种 植 区 , 小 区 随机 排列 , 共 设 9 个 处 到 


抽 干 时 气象 与 生态 环境 试验 基地 (35°33 


黄土 高 原 西部 丘陵 区 和 半 干 旱地 区 , 具有 典型 的 温带 大 陆 性 
季风 气候 L199。1971 一 2014 年 平均 降雨 量 388 mm, 冬季 和 夏季 月 平均 降雨 量 分 别 为 20~80 mm 和 150-270 


年 降雨 量 的 86.9%,， 蒸发 强烈 ,年 均 潜 在 蒸发 量 达到 1500 


cm 土壤 平均 容重 为 1.38 g.cm3， 田 间 持 水 量 为 25.6%, 饱和 
含水 量 为 43.87%, 永久 凋 芋 系数 为 6.7%, 地 下 水 埋 深 10.4 m, 土壤 水 与 地 下 水 的 水 力 联系 微弱 。 当 地 耕 
作 制 度 为 1 年 1 XA, 主要 种 植 作物 有 春小麦 (Triticum aestivum)、 燕 麦 (Avena sativa) fl H48 È (Solanum 
tuberosum) 等 ， 主要 种 植 牧 革 有 紫花 首 和 菠 (Medicago sativa) fl £L S 5t (Onobrychis viciaefolia) S . 

试验 采用 田间 垄沟 集 雨 覆盖 种 植 技 术 ， 以 燕麦 (4vena sativa) 为 指示 作物 ， 垄 覆盖 作为 集 十 区 , 沟 无 履 
HE 3 次 。3 种 覆盖 材料 分 别 


R(3 种 沟 垄 比 x3 RIS MED, 8 


为 生物 可 降解 膜 、 普 通 塑 料 膜 和 士 壤 结 皮 , 3 种 沟 克 比分 别 为 60 cm : 30 cm、60 em : 45 cm 和 60 cm : 60 


cm( 沟 宽 : 垄 宽 )。 生 物 可 降解 地 膜 和 普通 塑料 膜 厚 度 均 为 0.08 mm; 了 


区 为 人 工 原 土 夯 实 ,， 经 风 吹 雨 打 形 成 


自然 土壤 结 皮 。 土 蕉 、 生 物 可 降解 地 膜 垄 和 普通 地 膜 垄 的 代表 符号 分 别 为 SR、BMR 和 CMR。 根据 当 地 
长 为 10m, 每 1 小 区 有 4 条 垄 和 3 条 沟 。 试 验 设计 示意 图 如 


种 植 经 验 , 集 雨 垄 的 坡度 为 40"， 高 为 25 cm, 
图 1 所 示 。 相 关 种 植 管理 方法 按 前 期 


统 土 壤 水 分 时 空 动态 影响 不 显著 0619， 故 本 文 结果 和 讨论 部 分 月 


种 试验 沟 歼 比 的 平均 值 。 
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究 者 9 所 述 进行 。 由 于 本 试验 的 3 PREY HAZRET A 
到 的 试验 数据 , 均 为 同 1 覆盖 材料 下 ,3 


土著/ 普通 膜 欧 / 生 物 可 降解 腊八 | 
|. Ridge with compacted soil 1 
| or covered with common plastic film | 
or biodegradable mulch film. 1 


RRR 降雨 Phting z 
Evapotranspiration Rainfall L 60 cm or 45 Jd a Omosa E 
' , 1. 燕麦 种 植 示意 图 
Plant interception Fig 1.1 Schematic diagram showing ridge-furrow rainwater harvesting for oats production 
pi 
植物 截留 表 1.1 低 沟 集 雨 种 植 燕麦 试验 设计 
Table 1.1 Experimental design for oats production in ridge-furrow rainwater harvesting 
处 理 i925 c (em) zig fm?) 沟 面积 m2?) 小 区 面积 (m”) 玖 覆盖 方式 
| Treatment Furrow : Ridge Ridge area Furrow area Plot area Ridge mulch style 
SR30 60:30 13.5 18 31.5 土壤 结 皮 
SR45 60:45 20.25 18 38.25 Soil in 
i oil crust 
MAN H i SR60 60:60 27 18 45 
j SA i 
d Y) A) zi ! 3 B .5 31.5 
3 X 根系 层 深 Dr， BMR30 60:30 13 18 1 可 降解 地 膜 覆盖 
Pp \ Root depth i BMR45 60:45 20.25 18 38.25 * ] 
RIR Eh SUE s Ro a , T: as Biodegradable mulch flm | 
Effective porosity : CMR30 60:30 13.5 18 31.5 ETTTLIS | 
VERMES CMR45 60:45 2025 18 3825 — 
IR ES 漏 Leakage CAES ED us B FE Common plastic film 


图 1 垄沟 集 雨 系统 土壤 水 分 动态 过 程 及 其 随机 模拟 试验 设计 示意 图 
Fig. 1 Schematic diagrams showing experimental design for stochastic simulation of soil moisture dynamics in ridge-furrow 
rainwater harvesting 
试验 期 降雨 量 数据 由 试验 基地 自 记 雨量 计 测 定 。 在 燕 考 播种 前 (4 月 10 HEA 收获 后 (8 月 20 HÆ 


BH) SCS-CN 模型 反 推 


右 ) 和 降雨 后 (降雨 量 >5 mm) 测 定 土 壤 含 水 量 , 土壤 含水 量 采 用 烘 干 法 (105 C, 10 WE, 测定 深度 为 140 
cm， 分 层 深 度 为 20 em, KUR 7 个 土 层 的 土壤 含水 量 ， 
REKER 3 个 样 点 平均 值 。 集 雨 歼 径流 量 由 降雨 量 数据 基于 美国 水 土 保持 局 研 


1 小 区 随机 在 沟 内 选取 3 个 样 点 , 同 1 层次 土 


确定 20。 土 壤 容 重 采用 环 刀 法 测定 , 测定 深度 140 cm， 分 层 深 度 为 20 cm, 每 层 3 个 重复 ， 取 均值 。 根 系 
层 深度 通过 实地 调查 燕麦 根系 生物 量 分 布 范 围 测 定 趾 。 同 期 气象 数据 由 邻近 的 试验 基地 气象 观测 站 获得 。 
1.2  Laio 土壤 水 分 动态 随机 模型 
土壤 水 分 随机 模型 的 理论 基础 是 物质 平衡 原理 : 单位 时 间 内 土壤 含水 量 的 变化 等 于 土壤 水 分 输入 项 
和 水 分 损失 项 的 甜 。 基 于 前 期 研究 者 2 对 土壤 水 分 随机 模型 的 研究 ，Laio 等 0 通过 引进 两 个 临界 土壤 含 
水 量 (土壤 吸湿 系数 和 土壤 凋 蓉 系 数 ), 在 空间 一 点 上 建立 了 时 间 尺 度 为 1d 的 土壤 水 分 平衡 方程 ,具体 模 
型 (Laio 模型 ) 可 表述 为 : 
ds(t) 
t dr 
式 中 : on IERIURE;, s 为 土壤 饱和 度 ; z, 为 根系 层 深 度 (cm); s() t 时 刻 的 土壤 相对 湿度 
(0xs(x1), 即 s()=00)/n， HPOO 为 上 时 刻 的 土壤 体积 含水 量 (9%); sH: ze PEERS N38 AE Ccm d^), 
表示 实际 到 达 土 壤 的 降雨 量 部 分 ， B ops(t);t] 2 RG)—1() - Qls(0;0], APRO, I) 和 Cs(D; Mie 
水 速率 (cmd0D) 、 植 物 截留 率 (cmd0D0 和 地 表 径 流 率 (cm'd0D; AiR M A Js TJ pede I, B 
X[s(D1= EOIs, 其 中 EBO 792880 s SE Cmd, MO 为 深层 渗 漏 率 (cm.dD)。 
pe 34 EE RS BEC OE RR Fe] E RKB JT BE] ERKI R UU £r Je ERKA B D LC Se v B5 mu 
p P2122。 由 于 降雨 是 随机 过 程 ， 故 需 建立 土壤 水 分 概率 密度 函数 (Soil Moisture Probability Density Function) 
求解 土壤 水 平衡 过 程 ( 式 1)。 通 过 将 Laio 模型 的 各 土壤 水 分 损失 过 程 转化 为 查 普 曼 - 柯 尔 莫 哥 洛 夫 前 进 方 
程 (Chapman-Kolmogorov Forward Function) 可 分 析 求 解 导 出 土壤 水 分 概率 密度 函数 ， 其 具体 表达 式 
7.1121-22] 为 : 


nz = q[s(t):r]— ZSO] = RO - UG) + lst); n] - ELO + HAsO] (1) 


T 


lo PIE | 
Coors| 575. | " | s, <S <SS 
Ny Sw =S, 
A (s' -s, ) 
5 Hh 
MM am ， s, <SSS (2) 
Ny, 1y S 一 Sw 
p(s) = guias 
C 3e (eoa) 7 Thy * 
Le n E ; $ «s&€ Ske 
7 7, 
49) À' yi A (s -s Pt jm 
€ [7G B)s* Bs, FA n -5)y ps PT- 由 "^ 
= £5 d — UE m Sk <s< 1 
, 7 (n - m)e^* - me^ 7), 
SE; AME 1 ` p He I 
式 中 : 常数 c 由 | pods =1 求 解 得 出 ,其 他 模型 参数 含义 及 取 值 范围 见 表 2 所 示 , 各 处 理 待 优 参数 值 由 蒙 


特 卡 洛 试 错 法 确定 。 模 型 相关 假设 条 件 及 详细 推导 过 程 详 见 文献 0。 
表 1 垄沟 集 雨 系统 Laio 模型 参数 取 值 范围 及 其 不 同 处 理 下 的 待 优化 值 


Table 1 Values ranges and initial values of parameters in Laio model in ridge-furrow rainwater harvesting system 


随机 模型 参数 "m 参数 取 值 范围 不 同 处 理 参数 待 优化 值 
g "YS 
属性 Parameters of Laio model 单位 Values ranges of parameters Initial values of parameters 
Property p "n Units 
fj * Symbol 意义 Meaning UBV LBV SR CMR BMR 
土壤 了 筷 阶 度 
n ERARE mm? 0 1 0.58 0.58 0.58 
Soil porosity 
土壤 孔隙 大 小 分 布 参 类 
B RURANA 3 / 10.00 20.00 14.8 14.8 14.8 
Parameter describing drainage in soil 
土壤 活动 层 深度 
Z, WERE cm 140 140 140 140 140 
Depth of rooting zone 
土壤 饱和 导 水 率 
Ks RAA TKE cm. d! 10.00 30.00 19.5 19.5 19.5 
Saturated hydraulic conductivity of soil 

土壤 Bağ 
p 吸湿 系 
Soil Sh CIEN mim? 0 0.10 0.026 0.026 0.026 


Absorption coefficient of soil moisture 


DEEST - 


Sw m? 0.10 0.40 0.25 0.25 0.25 
Degree of soil saturation at wilting point 
7] Z ME 6H 
s TO REOR mm? 0.40 0.70 0.56 0.56 0.56 
Soil saturation at stress point 
间 持 水 率 
Sfo 司 持 水 率 m? 0.70 1 0.72 0.72 0.72 
Degree of soil saturation at the field capacity 
凋 萎 系数 对 应 土壤 蒸发 
Ey 辣 蔡 系 数 对 应 土壤 燕 发 cm 0 0.02 0.012 0.012 0.012 
Estimated rate of evaporation from soil 
FEL ER EA N 
作物 A f MRE RR . cm 0 0.50 0.128 0.130 0.131 
Interception capacity of vegetation 
Crop E3 J-& ir. E 
Emax ] ROGER C cm 0 1.00 0.610 0.570 0.560 
Maximum evapotranspiration 
生长 季 平 均 降 水 量 同 初 始 值 
ES a cmd! 0.5949 0.8166 0.7673 
降雨 Mean rainfall depth of growing season Same as the initial value 
Rainfall 生长 季 平 均 降 水 频率 
À 生长 季 平 均 降 水 频率 di 0.33 0.33 0.33 0.33 0.33 


Mean rainfall frequency of growing season 


表 中 模型 参数 取 值 范围 由 相关 实测 数据 和 文献 资料 "19 获 得 ; UBV: 参数 取 值 下 限 ; LBV: 参数 取 值 上 限 ; SR. CMR 和 BMR 分 别 表示 土 垄 、 普 
通 膜 共和 生物 可 降解 膜 垄 。 下 同 。The range of model parameters in the table is obtained from the relevant measured data and literature (e.g., [11,16] ); 
UBV: upper bound value of model parameters; LBV: lower bound value of model parameters; SR, BMR and CMR were ridges with compacted soil, covered 
with biodegradable mulch film and common plastic film, respectively .The same below. 


13 ”模型 参数 敏感 性 分 析 原 理 

模型 参数 的 敏感 性 分 析 是 研究 参数 变化 所 引起 的 模型 响应 ,是 模型 参数 不 确定 分 析 的 重要 内 容 之 一 
也 是 研发 和 评价 模型 不 可 缺少 的 重要 环节 中 。 同 时 ,参数 敏感 性 分 析 有 助 于 深入 理解 模型 的 特性 并 改进 
模型 结构 的 稳定 性 多。 为 不 失 一 般 性 , 将 土壤 水 分 概率 密度 函数 表示 为 : 

p=f(%,%,L Apta) (3) 
式 中 : p 为 土壤 湿度 概率 密度 ; x 为 第 ; 个 影响 因素 ; n 为 影响 p 的 因素 个 数 。 
当 所 有 因素 都 发 生变 化 , IIH, x. L, x 79x, x, L, x, 1G Ax, AX, L, 
Ax Wf, xx Ax, xox-Ae, L, x-x-tAx, Ml Sie SERE WS BE p 也 发 生 相 应 变化 , BH p AED p', 
则 可 用 Ap = p'— p 表示 所 有 因素 变化 共同 造成 p 的 变化 量 。 利 用 多 元 函数 的 泰勒 展开 式 为 : 

Ang OP ge EP as +L A (4) 

Ox, Ox, Ox; 

式 中 : Op/0x, 73 x, 对 pp 的 偏 导 数 ，Ax 为 x 的 变化 量 。 如 果 只 有 x 因素 改变 , 其它 因素 都 不 发 生变 化 ， 即 
Ax 0 ，A% =0，1zi， 则 土壤 湿度 概率 密度 p ZEWN Ap, RE A, X p 的 影响 值 ， 可 表示 为 
Ap, ~( 印 /x) Ax, 显然 p 的 综合 变化 量 为 Ap=》"Ap, o 

定义 的 变化 量 与 p 因素 zx 的 变化 量 之 比 为 p 对 x 的 敏感 度 A4， 则 A=(Ap,/p)/(Ax/x,)， 其 表示 第 i 个 
因素 二 变化 1 个 百分点 , 将 引起 土壤 湿度 概率 密度 pp 改变 A 个 百分点 。4 为 正 表 示 Ap, 与 Ax 的 变化 方向 
相同 。 JA RK, 表明 因素 x 对 土壤 湿度 概率 密度 p 的 影响 越 大 ， 即 p 对 越 敏 感 ， 可 认为 x 为 敏感 因素 。 
敏感 度 4 计算 公式 如 下 : 

A M/k [EAk dp x 

| Ax, /x AX Lg, Ox, k 

在 获取 垄沟 集 雨 系统 Laio 模型 各 参数 敏感 度 时 ， 首 先 借 助 Mathematica 8.0.4 软件 对 模型 各 参数 敏感 
E A 的 通 式 进行 求 取 ,然后 选 定 土壤 含水 率 * 的 范围 (0<s<1), 采用 表 1 中 不 同 处 理 的 待 优化 参数 值 ， 
计算 各 处 理 对 应 参数 的 敏感 度 ， 确 定 不 同 处 理 下 模型 敏感 性 参数 ， 同 时 基于 SOM 神经 网 络 聚 类 法 对 垄沟 
集 雨 系统 Laio 模型 参数 敏感 性 进行 分 类 ， 并 探讨 同一 敏感 性 类 别 的 参数 特征 及 不 同 敏感 性 类 别 形成 的 主 
要 原因 。 


14 ”模型 参数 优化 算法 原理 


1.4.1 改进 单纯 形 法 寻 优 原理 

单纯 形 法 (Simplex Method, SM) 也 称 可 变 多 面体 搜索 法 , 是 1 种 局 部 搜索 方法 , 由 Nelder 和 Mead 于 
1965 年 提出 RI。 其 基本 思想 是 : 在 NN 维 欧 氏 空间 中 取 N+1 个 点 构成 初始 单纯 形 , 假设 每 个 点 的 坐标 矢量 
H XX L, Xy EEX Nel 个 点 处 目标 函数 值 的 大 小 , 求 出 函数 值 最 大 点 X, 对 剩余 点 
X,X,L,X, X, oL Xy JE X, 的 映射 和 AEH x, STAR xX, ,构成 新 的 单纯 形 ， 反 复 迭 代 ， 使 
其 顶 点 处 的 函数 值 逐步 下 降 从 而 逼近 目标 函 数 ， 即 采用 扩张 、 收 缩 、 压缩 和 扬弃 这 4 种 方法 对 各 点 进行 
处 理 ,直至 满足 收敛 性 指标 , 搜索 到 最 佳 点 。 

本 文 的 改进 单纯 形 法 (Improved Simplex Method, ISM) 是 在 基本 单纯 形 法 的 基础 上 对 步 长 作 适 当 修 改 
得 到 的 寻 优 方法 R99。 设 定 最 佳 点 (函数 值 F 最 小 )、 次 佳 点 (函数 值 F 次 小 ) 和 最 差点 (函数 值 F 最 大 ) 的 坐标 
矢量 xX, X, x. 的 函数 值 分 别 为 已 ,已 和 已 ,映射 点 和 的 函数 值 为 忆 ， 于 是 根据 反射 值 己 可 作 不 同 
的 选择 : 

(DD) 如 果 巨 < 所 ,说 明 反 射 方向 正确 , 求 扩张 点 坐标 矢量 X =X +42x(X,_X,) ， 其 中 
X.=(X+X E +X, +% L+ f FECE, WI SKKE, Hi x, REX, ; 否则 扩张 失败 , 保 
Hx. 

(2) 如 果 石 << 五 ， 则 不 扩张 也 不 收缩 , 保留 Xx， ; 

(G) 如 果 已 , < 已 ， 说 明 反 射 方向 不 正确 ， 需 进行 收缩 ， 求 收缩 点 坐标 矢量 和 = X, -05*(xX,-xX,).) 
X, 代替 和: 

(DWR F, <F.<F,, 求 收缩 点 坐标 矢量 X, =X, +0.5*(X,-X,); HX, RË X,- 

按 上 述 x, 的 计算 结果 ,从 新 构造 新 的 单纯 形 ， 作 收敛 性 检验 ， 反复 迭代 直到 满足 终止 条 件 为 上 ,ISM 
迭代 是 否 结束 ， 需 根据 事先 给 定 的 条 件 来 判断 通常， 可 用 系统 啊 应 标准 偏差 S. 与 系统 响应 平均 值 到 之 比 
COV (偏差 系数 ) 作 为 判 据 , X? COV <0.5% ， 可 终止 优化 。 

1.4.2 ”粒子 群 优化 算法 寻 优 原 理 

粒子 群 优化 算法 (Particle Swarm Optimisation，PSO) 是 1 种 基于 群 智 能 方法 的 演化 计算 技术 ,由 
Kennedy 和 Eberhart 于 1995 年 提出 B27 在 PSO 算法 中 , 每 个 优化 问题 的 潜在 解 均 可 看 作 是 搜索 空间 中 的 
1 个 粒子 。 每 个 粒子 均 有 1 个 被 优化 函数 决定 的 适应 值 (Eitness Value, FV), 同时 有 1 个 速度 决定 的 移动 方 
癌 和 距离 ， 可 通过 跟踪 当前 最 优 粒 子 的 个 体 极 值 和 整个 种 群 的 全 局 极 值 迭 代 搜 索 最 优 解 。 

PSO 算法 的 基本 数学 描述 为 P3: 设 在 1 个 DD 维 的 目标 搜索 空间 中 ,， 用 个 粒子 组 成 上 个 群体 ,其 中 
第 ;个 粒子 表示 为 1 个 万 维 向 量 X=(xuxuL ,zj(i=1L2L N), 每 个 粒子 的 位 置 就 是 1 个 潜在 的 解 。 将 
X, 代 入 1 个 目标 函数 均 可 求 出 其 FV, 根据 FV 的 大 小 衡量 解 的 优 务 性 ,第 i 个 粒子 的 移动 速度 也 是 D 维 
向 量 ， 记 为 Y=(woviL wo)(i=42L ,N) 。 记 第 i 个 粒子 迄今 为 止 搜索 到 的 最 优 位 置 为 
P-(py. PaL ,po)(i=12L ,N)， 整 个 粒子 群 迄今 为 止 搜索 到 的 最 优 位 置 为 忆 =(P PsL Py) 。 则 可 | 
下 述 公式 对 粒子 的 移动 速度 和 位 置 进行 更 新 : 

yi +1)=v; (t)w+ nc [py —x 一 入 (t) «nc, | ei =% (0) (6) 

X; (r1) =X; (t)+v; (t+1) (7) 

AP, rn 和 是 服从 U(0,1) 分 布 的 随机 数 ; c RE c, 是 随机 因子 , 通常 取 a=c=2; w 是 惯性 因子 。 在 
每 1 维 中 ,粒子 均 有 1 个 最 大 限制 速度 W.， 如 果 某 1 维 的 速度 超过 V,, ,那么 这 1 维 的 速度 就 限定 为 Vs,， 
HI [v] V... 

143 ”混合 粒子 群 优 算法 寻 优 原理 

为 提高 优化 算法 在 全 局 和 局 部 意义 下 的 搜索 能 力 和 收敛 效率 ， 以 粒子 群 优 算法 流程 为 基础 ， 引 入 改 
进 单纯 形 搜索 方法 构成 混合 粒子 群 优 算法 (Hybrid Particle Swarm Optimization, HPSO)29]， 即 在 每 1 次 迭代 
中 先 用 PSO 算法 对 群体 进行 全 局 寻 优 ,然后 通过 ISM 算法 对 粒子 群 中 部 分 精英 粒子 在 其 较 优 解 领域 内 进 
行 局 部 搜索 ,寻找 更 优 解 。 
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HPSO 算法 的 具体 流程 9 为 : 对 每 次 迭代 中 经 过 PSO 优化 后 的 所 有 P 个 群体 粒子 (P 为 种 群 数 ) 按 适 


应 度 优 劣 值 进行 排序 , 选 出 适 
点 中 选 出 最 优 响应 的 9-1 个 顶点 马 , 忆 工 X，， 
“伸缩 ”映射 生成 XX' 顶点 ,构造 


出 新 的 “单纯 形 ” 


每 1 新 粒子 适应 度 值 ， 从 ， 


子 群 进化 计生 


选 出 最 优 响应 的 粒子 蔡 换 原先 群体 中 的 最 伯 
成 下 1 代 新 群体 。 接 着 比较 新 个 体 适 应 度 值 与 其 自 


并 计划 


达到 所 人 允许 的 最 大 代数 7 或 评价 值 小 于 


RM mee 
给 定 精度 =， 则 输 晶 


RTS, MEM SEREK, S/P 值 在 109-20% 之 间 为 宣 。 


1.5 


Laio 模型 转换 与 参数 优化 


应 度 值 最 优 的 前 8 个 粒子 构造 1 个 包含 8 个 顶点 的 “单纯 形 ” 图 形 ， MS 个 顶 
5 一 1 个 顶点 形 心 X， 由 剩余 顶点 X, 通过 形 心 X 

图 按照 上 法 多 次 重复 产生 S 个 新 粒子 。 然 后 计算 更 新 后 
个 体 ， 并 与 原 群 体 中 剩余 个 体 构 
最 佳 pbest 位 置 。 若 粒 
最 优 结果 pbest, ERAR. — 


非 线性 模型 参数 优化 问题 是 在 模型 结构 (表达 式 为 y= (6) e, e ~ N (0,5?)， 其 中 ，x 为 系统 输入 ， 
) 已 确定 的 情况 下 ,根据 已 知 


y 为 系统 输出 ，9 为 模型 待 优化 参数 集 ，e 是 均值 为 0、 方差 为 2 的 白 噪声 。 
的 输入 输出 数据 (x,y,)(i=1,2,L ,n)， 


求 出 存在 于 非 线性 模型 中 的 参数 集 9 


， 即 求解 偏差 平方 和 


F(0) =F (x-f (x0) 为 最 小 的 6 值 59。 为 不 失 一 般 性 ,对 垄沟 集 雨 系统 Laio 模型 进行 参数 优化 时 ， 


根据 已 知 数据 (s.s p(s),)(i=1,2,L n). X H 
f£ o 即 为 模型 参数 优化 值 。F (9) 越 小 , 表明 优化 参数 集 
优化 求解 需 综合 考虑 模型 参数 的 限 币 
用 于 非 线性 、 不 可 微 和 多 峰值 的 复杂 函数 优化 问题 的 求解 6 
分 别 对 Laio 模型 


算法 PSO 和 混合 搜索 全 


法 HPSO 4 


| 条 件 ， 故 传统 全 


E JB FEAR; I 7 > 


法 很 难 对 其 进行 求解 。 


ISM, PSO 和 HPSO 算法 收敛 效率 及 优化 参数 有 效 性 进行 分 析 和 验证 。 


1.6 


模型 的 验证 与 评价 方法 
通过 比较 研究 区 不 同人 处 型 


Measure) 指 数 B23， 其 表示 两 目 
— 24 > 


A t A 
式 中 ，4 表示 目标 | 


ERI 


1 线 工 下 方 的 面积 
A2 下方 的 公共 面积 ( 见 图 2) SA, 目标 | 
Z, 一 致 性 越 差 ，A4,, 越 小 , CM 指数 越 小 。CM 指数 
目标 曲线 2 一致 性 最 佳 ; 当 CM=0 Hf, 表明 目 


实测 土壤 水 分 概率 密度 函数 曲线 与 ISM, 
拟 曲线 的 特征 , 评价 以 上 3 种 算法 对 Laio 模型 参数 优化 的 有 效 性 。 验 说 
| 线 的 一 致 性 程度 ,计算 方法 如 下 : 


| 线 1 与 目标 | 


y 


， 册 表示 目标 曲线 2 下方 的 面积 


(a) 


AX 


明 越 不 可 靠 。 


标 函 数 居 (0O)=min> (p( (s), - f (s. 8) 进行 求解 ， 所 得 参数 
目标 函数 天 (b) 的 


由 于 众多 现代 非 传统 算法 可 
HH。 故 本 文选 择 局 部 搜索 算法 ISM、 全 局 搜索 


(8) 


4 ,表示 目标 


曲线 1 和 目标 曲 
12k 2 至 性 越 好 ，4 RK, CM 指数 越 大 ; 反 
的 变化 范围 为 [0, 1], 24 CM=1 时 , 表明 目标 


偏差 平方 和 (9) 的 最 小 值 进行 优化 求解 ， 并 分 别 对 


PSO 和 HPSO 算法 优化 参数 值 模 
E 评 价 指标 选择 CM(Consistency 


H 线 1 和 


(b) 


图 2 一 致 性 测度 (CM) 指 数 定义 示意 图 


Fig.2 Schematic diagram showing definition of consistency measure(CM) 


标 曲 线 1 和 目标 曲线 2 之 间 不 存在 一 致 性 ( 见 如 2(b))。 
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2 结果 与 分 析 
2.1 Laio 模型 参数 敏感 性 分 析 
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3 ”垄沟 集 雨 系统 Laio 模型 各 参数 敏感 性 分 析 结 果 
Fig.3 The sensitivity degree of parameters in Laio model in ridge-furrow rainwater harvesting system 
采用 表 1 中 不 同 处 理 待 优化 参数 值 ， 借助 MATLAB R2010a 软件 ， 计 算得 出 不 同 土壤 含水 率 下 
(S,€«s€s,. s,«s«s . s «scs, Ts, «s«l)Laio 模型 各 参数 的 敏感 度 (图 3)。 从 图 3 可 以 看 出 , 在 不 


同 土壤 含水 率 和 不 同 处 理 条 件 下 ,模型 参数 敏感 性 表现 不 同 。 总 体 而 言 ， 6 Ko Asp 对 模型 的 输 


影响 甚 微 A 和 Ema 对 模型 的 输出 影 


啊 较 小 ; e RI s, 对 模型 的 输出 


LEVA 
R 


DG 


响 较 显著 ,为 影响 垄沟 集 雨 系统 土 


壤 水 分 动态 随机 模拟 的 最 重要 因子 ; 土壤 含水 率 为 5, <s<s, 时 ，n 、Z,、E, 和 4 对 模型 的 输出 影响 较 显 


著 ，* 对 模型 的 输出 影响 甚 微 ,而 在 其 


他 土壤 含水 率 条 件 下 ， 其 均 对 模型 的 输出 影响 较 小 ; 各 处 理 对 应 


模型 参数 敏感 度 随 士 壤 含水 率 的 变化 趋势 基本 一 致 。 在 相同 土壤 含水 率 条 件 下 ， 处 理 SR 对 应 模型 敏感 性 
参数 的 敏感 度 显 著 大 于 处 理 CMR 和 BMR, 处 理 CMR 与 BMR 之 间 差 异 不 显著 。 在 相同 处 理 条 件 下 , 低 


土壤 含水 率 对 应 模型 敏感 性 参数 的 敏感 度 显 著 大 于 高 含水 率 条 伯 


F 下 的 敏感 度 。 同时 ， 除 参数 、w 和 4 对 


应 敏感 度 为 正 值 外 ， 其余 各 模型 敏感 性 参数 的 敏感 度 均 为 负 值 ， 且 均 为 土壤 特征 参数 。 
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图 4 基于 SOM t kaai Laio 模型 参数 不 同 分 类 结果 


Fig.4 Different classification results of parameters in Laio model based on SOM neural network clustering method 


采用 SOM 神经 网 络 聚 类 法 对 垄沟 集 雨 系统 Laio 模型 各 参数 敏感 性 进行 分 类 (图 4)。 从 图 4 可 以 看 出 ， 


将 13 个 模型 参数 分 为 3 类 时 , 不 同 土壤 含水 率 条 件 下 的 分 类 结果 存在 相同 之 处 , 即 6 K, 和 si; 均 为 弱 敏 
感 参数 ; % 和 4 均 为 中 等 敏感 参数 ，w M s, 均 为 强 敏感 参数 。 同 时 ,不 同 土壤 含水 率 条 件 下 w% > Ar Emo 
和 sw, 的 分 类 结果 存在 差异 。 通 过 综合 考虑 不 同 土壤 含水 率 条 件 下 模型 各 参数 敏感 性 分 类 结果 , 将 垄沟 


集 雨 系统 模型 参数 敏感 性 可 分 为 3 类 : 第 一 类 为 弱 敏 感 参数 , 包括 参数 6 Ko Ms; 第 二 类 为 中 等 


PRSA, UBA, Ema ns Zo E AAMs, 第 三 类 为 强 敏 感 参数 , 包括 参数 x 和 s,。 这 与 姚 淑 
起 等 四 的 分 类 结果 基本 一 致 ，x 和 $s, 均 为 强 敏 感性 参数 ，K, 为 弱 敏 感性 参数 ; 而 与 任 庆 福 0 的 分 类 结果 


不 太一 致 ， 


,在 任 庆 福 分 类 中 为 弱 敏 感性 参数 , 但 在 本 文 分 类 中 为 强 敏 感性 参数 ,造成 这 种 差异 的 原因 


可 能 EARN THLR A Ke Fh deo 行 综 合 考虑 。 同 时 , 不 同 研究 区 土壤 质地 、 植 被 和 
气候 特征 之 间 存 在 的 差异 ,也 可 能 是 造成 模型 参数 敏感 不 一 致 的 原因 ， 故 应 根据 研究 区 的 具体 环境 条 件 
对 模型 参数 敏感 性 做 相应 分 析 。 


100 


模型 模拟 结果 的 


的 增 大 引起 p(s) 的 非 线性 减少 ; 参数 w 在 0.135~1 cm:d-! 范 围 内 变化 时 ，w 的 增 大 引起 p(s) 先 非 线性 增加 


后 非 线性 减 小 , Ha 73 0. 


Qi 
ps) 


图 5 参数 a 和 5,, 对 垄沟 集 雨 系统 土壤 水 分 概率 密度 函数 p(s) 的 影响 


Fig.5 The influence of parameter e and parameter s, on p(s) in ridge-furrow rainwater harvesting system 


由 于 强 敏感 性 参数 wx F s, 对 模型 模拟 结果 的 影响 显著 高 于 其 他 参数 ， 故 需要 进一步 探讨 xc 和 5s, 对 


影响 。 以 处 理 CMR 为 例 ， 当 其 他 参数 不 变 时 , 讨论 < RI s, e ELI] CHAR LOG] 2E] E Rf 2 
统 土壤 水 分 概率 密度 函数 p(s) 的 影响 (图 5)。 由 图 5 可 知 ， 参 数 w 在 0.026-0.135 cmd! 范围 内 变化 时 ，cw 


Wi, 参数 的 变化 能 


192 cmd! If, p(s) 达到 最 大 值 。 同 时 ,土壤 含水 率 在 0.026~0.25 cm?-cm? 范围 
起 p(s) 显著 波动 ,说明 土壤 水 分 概率 密度 函数 p(s) 对 的 敏感 性 在 低 含水 率 下 


更 明显 。 表 明 在 干旱 半 干 旱 区 ,垄沟 覆 膜 集 雨 系统 可 将 无 效 降水 (小 于 5 mm) 转 化 为 有 效 水 分 贮存 于 土壤 ， 
显著 提高 土壤 含水 率 0331。 人 参数 , 的 变化 范围 为 0.026~0.56, 在 土壤 含水 率 区 间 内 ，p(s) 呈 非 线性 增加 趋 
A, 且 土 壤 含水 率 越 大 ，y*, 变化 引起 的 p(s) 变化 越 大 , 说 明 土 壤 水 分 概率 密度 函数 p(s) 对 s, 的 敏感 性 在 
高 含水 率 下 更 明显 ,表明 垄沟 履 膜 集 雨 系统 能 在 高 效 雨 水 收集 的 基础 上 ,使 作物 摆脱 干旱 胁迫 ,显著 提高 
土壤 水 分 利用 效率 034。 
2.2 基于 ISM、PSO 和 HPSO 算法 的 模型 参数 优化 

根据 参数 敏感 性 分 析 结 果 及 参数 优化 原则 , 分 别 利 用 ISM、PSO 和 HPSO 算法 对 Laio 模型 中 的 n、 
Bo Ko o Spo Sh Sye En A, Em 11 个 参数 进行 优化 。 各 待 优 参 数 初 值 及 取 值 范围 见 表 1。 
在 参数 优化 过 程 中 , 各 算法 参数 设置 如 下 : 种 群 数 均 设 为 80; 最 大 迭代 次 数 为 1000; 适应 度 函 数 了 采用 
偏差 平方 和 形式 ( 见 1.5 中 ); ISM 算法 扩张 系数 y=2， 压缩 系数 6=05; PSO 和 HPSO 算法 速度 更 新 式 中 惯 
性 权 系 数 w 值 线性 递减 , 其 中 mw，=0.9，mwi 204; 学 习 因 子 取 为 c =c 202; HPSO 算法 中 用 于 改进 单 


纯 形 局 部 搜索 的 优秀 个 体 数 $=6, 扩张 系数 yx=2,， 压缩 系数 6=05; 整体 算法 每 进化 一 代 , 用 改进 单纯 
形 法 对 选中 的 优秀 个 体 搜索 100 步 ; 参数 Z,、w 和 14 均 根据 试验 实测 值 确 定 。 

为 确保 算法 的 优化 精度 和 有 效 性 , 在 算法 运行 的 过 程 中 符合 下 列 3 个 终止 准则 之 一 时 停止 计算 B51。 
准则 1: 两 次 迭代 的 最 优 目标 函数 值 对 应 的 参数 距离 小 于 给 定 精度 a «107 ; 准则 2: 两 次 迭代 的 目标 函数 
值 之 差 小 于 给 定 精度 e «107 ; 准则 3: 循环 的 最 大 次 数 已 达到 。 同 时, 为 尽 可 能 消除 算法 随机 性 对 各 算法 
比较 的 影响 ,每 个 优化 算法 均 随机 独立 运行 10 次 后 取 各 参数 平均 最 优 解 和 计算 效率 。 算法 运行 结果 和 优 
化 参数 验证 见 表 2 及 图 6 所 示 。 


表 2 基于 ISM、PSO 和 HPSO 算法 的 垄沟 集 雨 系统 Laio 模型 参数 优化 结果 
Table 2 Optimization results of parameters in Laio model in RFRHS based on ISM、 PSO and HPSO algorithm 


— 模型 参数 单位 ISM 算法 参数 优化 值 PSO 算法 参数 优化 值 HPSO 算法 参数 优化 值 
s Values after optimization using ISM Values after optimization using PSO Values after optimization using HPSO 
CC Parameters Units 
A | SR CMR BMR SR CMR BMR SR CMR BMR 
ph n m? m? 0.5696 0.5814 0.5811 0.5637 0.5863 0.5809 0.5634 0.5866 0.5804 
p / 14.8003 14.8011 14.8008 14.2084 13.5025 14.0283 14.2012 13.5481 14.0315 
Z, cm 140 140 140 140 140 140 140 140 140 
E K; cmd! 19.5004 19.5012 19.5007 15.5426 17.3944 17.5720 15.5407 17.3452 17.5699 
Sh m?-m? 0.0258 0.0263 0.0263 0.0221 0.0288 0.0243 0.0261 0.0279 0.0238 
= Sw m?-m? 0.2495 0.2490 0.2494 0.1801 0.2066 0.2017 0.1302 0.1866 0.1817 
s* m? m? 0.5803 0.5703 0.5698 0.6713 0.5611 0.5325 0.6910 0.5712 0.5341 
Sf m? m? 0.7210 0.7210 0.7200 0.8206 0.7439 0.7454 0.8803 0.7468 0.7435 
E, cm 0.0090 0.0063 0.0100 0.0152 0.0178 0.0182 0.0132 0.0167 0.0181 
A cm 0.1280 0.1293 0.1310 0.1438 0.1119 0.1171 0.1386 0.1123 0.1164 
Emax cm 0.6111 0.6712 0.5598 0.6115 0.6009 0.5316 0.6619 0.6002 0.5311 
a cmd! 0.5949 0.8166 0.7673 0.5949 0.8166 0.7673 0.5949 0.8166 0.7673 


À d! 0.3300 0.3300 0.3300 0.3300 0.3300 0.3300 0.3300 0.3300 0.3300 
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SMC 表示 土壤 含水 量 ; NCF 表示 正 态 曲线 拟 合 ; Spso 表示 PSO 优化 参数 模拟 曲线 ; Sim 表示 ISM 优化 参数 模拟 曲线 ; SHpso 表 示 HPSO 优化 参数 


模拟 


线 。SMC represents soil moisture content; NCF represents normal curve fitting; Spso represents PSO optimization parameters simulated curves; Sism 


represents ISM optimization parameters simulated curves; SHpso represents HPSO optimization parameters simulated curves. 


图 6 基于 ISM、PSO 和 HPSO 算法 的 Laio 模型 参数 优化 效率 和 优化 参数 有 效 性 验证 结果 
Fig. 6 Optimization efficiency and validity verification of parameters in Laio model in RFRHS based on ISM, PSO and HPSO 


从 3 LE: 
法 优化 参数 


全 局 寻 优 能 力 的 算法 获取 一 


优 能 力 ， 优 化 参数 的 有 


组 较 


algorithm 


的 有 效 性 明显 不 及 另外 两 种 算法 ， 其 对 模型 参数 初 
优 初始 值 ， 其 优化 结果 将 易 陷 入 局 部 最 优 。 


陷 。HPSO 算法 


能 力 , 而 且 具 有 较 快 的 收敛 速度 ， 
29 对 HPSO 算法 的 仿真 实验 结 
的 优化 性 能 ， 是 求解 优化 问题 的 一 种 有 效 的 和 


A 


较 强 ， 


优化 参数 的 有 效 性 明显 优 于 ISM 算 ; 


[EE XE SETA 


PA 


能 较 显 著 地 克服 ISM 


法 ISM、PSO 和 HPSO) 对 Laio 模型 参数 优化 结果 ( 表 2) 和 优化 效率 ( 


台 值 依赖 性 较 


后 期 收敛 速度 相对 较 慢 ,存在 有 早熟 、 
和 PSO 算法 的 优化 解 ， 


y» 


法 和 了 PSO 算法 相互 


论 相 一 致 ， 说 明 其 有 利 


法 。 


2.3 ”模型 参数 有 效 性 验证 与 评价 


Y, 


增强 全 局 和 局 部 意义 


采用 试验 区 2012—2013 4E FH [8] 2874 48 Rl] 2 £05 Ab 
实测 数据 资料 基于 3 种 算法 ISM、PSO 和 HPSO) 对 Laio 模型 参数 优化 值 ， 
水 分 概率 密度 函数 在 曲线 形状 (峰值 、 峰 值 出 现 的 位 置 、90% 置 信 区 间 ) 和 CM 


3 种 外 


法 对 Laio 模型 参数 优化 的 有 效 性 ,选择 较 优 模型 参数 优化 值 及 参数 优选 算 ; 


PSO 算法 


图 6) 可 以 看 出 , ISM 算 
[I 果 不 事先 用 具有 较 强 
有 较 强 的 全 局 寻 
局 部 收敛 的 缺 
不 但 具有 较 强 的 全 局 搜索 
存在 的 缺陷 。 这 与 陈 俊 风 等 
下 优化 结果 的 可 靠 性 和 算法 


Ui 


强 
d. 


~ 


ECSR、CMR 和 BMR) 降 雨 、 径 流 、 土 壤 含水 率 等 


比较 模型 模拟 与 实测 的 土壤 
指数 之 间 的 匹配 程度 , 评价 
去 类 型 ,图 6 和 表 3 为 3 


IAR ZE 


种 算法 的 垄沟 集 雨 系统 Laio 模型 优化 参数 有 效 怕 


与 评价 结果 。 


EX ul 


表 3 基于 ISM、PSO 和 HPSO 算法 的 垄沟 集 雨 系统 Laio 模型 参数 优化 的 性 能 比较 
Table3 Performance comparison of parameter optimization in Laio model in RFRHS based on ISM, PSO and HPSO algorithm 


算法 SR 的 p(s) 曲 线 检验 指标 CMR 的 p(s) 曲 线 检验 指标 BMR 的 p(s) 曲 线 检验 指标 
Consistency test indicators Consistency test indicators Consistency test indicators 
Algorithm 
CPV PP CI95% CM CPV PP CI9596 CM CPV PP CI9596 CM 
ISM 0.0194 0.0333 -0.1493 0.519 0.0069 -0.0286 -0.0434 0.830 0.0111 0.0588 -0.0437 0.685 
PSO 0.0126 0.0333 -0.0874 0.723 0.0068 -0.0286 -0.0490 0.957 0.0093 0.0001 -0.0210 0.864 
HPSO 0.0001 0.0001 -0.0629 0.805 0.0061 -0.0286 -0.0241 0.994 -0.0024 0.0001 -0.0193 0.903 
CPV: 峰值 ; PP: 峰值 位 置 CI95%: 95% 置 信 区 间 ; CM: 一 致 性 指数 。 CPV: the curve peak value; PP: the position of the peak; CI9596: the 


confidence interval of 9596; CM: the consistency measure. 


从 图 6 和 表 3 可 以 看 H 
位 置 ， 描述 出 土壤 水 分 概率 密度 函数 的 


H, 3 种 算法 优化 模型 参数 值 均 能 比较 准确 地 描 


10%，CM 指数 均 大 于 0.5, 说 明 3 $05 
好 。 同 时 , HPSO $ 
PSO $ 


水 分 概率 密度 模拟 精度 和 效率 , 故 HPSO 算法 可 作为 垄沟 


法 优化 参数 值 对 垄沟 集 雨 系统 


zu 


[T 
Bj 


昌 线 的 形状 ， 撒 


和 捉 到 峰值 的 


ERRE, H CPV. PP 和 CI95% 实 测 值 与 模拟 值 相 对 误差 均 小 于 


上 : 壤 水 分 概率 密度 函数 模拟 效果 较 


法 优化 参数 的 模拟 效果 (CM 指数 均值 为 0.901) 优 于 PSO 算法 (CM 指数 均值 为 0.848), 
法 优 于 ISM 算法 (CM 指数 均值 为 0.678)， 日 HPSO 算法 优化 参数 的 收敛 速度 (进化 代数 均值 为 285) 
均 快 于 PSO 算 法 (进化 代数 均值 为 503), 说 明 HPSO 算法 优化 模型 参数 值 可 较 显 著 提高 垄沟 集 昌 


系统 
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雨 系统 土壤 水 分 动态 随机 模型 参数 优选 的 较 


优待 选 方案 。 但 是 ， 相关 研究 表明 9, 全 局 算法 与 局 部 算法 相 混 合 得 到 的 混合 优化 算法 ， 尽管 可 以 提高 后 
部 收敛 速度 和 性 能 ,但 也 加 剧 了 陷入 局 部 极 小 的 可 能 。 因 此 , 构建 能 显著 提高 局 部 搜索 能 力 ， 且 能 高 概率 
穆 索 全 局 最 优 解 的 混合 优化 算法 ， 有 待 进一步 深入 的 研究 。 
3 结论 

本 文采 用 多 因素 敏感 性 分 析 法 ,对 半 干 旱 区 芍 沟 集 雨 系统 Laio 土壤 水 分 动态 随机 模型 参数 的 敏感 性 


^ 


结论 如 下 : 


进行 了 分 析 和 分 类 。 同 时 , 基于 3 种 算法 
进行 了 优化 和 优选 ， 并 利用 实测 数据 资料 


| 对 3 种 个 


1) 通 过 综合 考虑 不 同 土壤 含水 率 条 们 


(ASM, PSO 和 HPSO), XJR n 
法 优化 参数 的 有 效 性 进行 了 验证 和 评价 。 得 到 的 主 


F 下 模型 各 参数 的 敏感 性 ,发 现 玲 


有 系统 Laio 模型 的 13 个 参数 


p: 
女 


4k 
ZR 


沟 集 雨 系统 Laio 模型 涉及 的 13 


个 参数 中 , 平均 年 降水 量 和 调 萎 系数 5 对 模型 输出 影响 最 显著 , 为 强 敏 感性 参数 ; 


2) 强 敏感 性 参数 a M s, 在 一 定 含水 率 


a 


e 


在 低 含水 率 下 更 明显 ， 而 对 参数 % 的 敏感 性 在 高 含水 率 下 更 明显 ; 


3) 三 种 入 


围 内 变化 时 , 土壤 水 分 概率 密度 函数 p(s) 对 参数 a 的 敏感 性 


法 的 优化 参数 值 均 能 较 好 模拟 垄沟 集 雨 系统 土壤 水 分 概率 密度 函数 ， 其 中 各 处 型 


PP 和 CI95% 实 测 值 与 模拟 值 相对 误差 均 小 于 10%, CM 指数 均 大 于 0.5; 


4)HPSO 算法 优化 参数 的 模拟 效果 和 收敛 速度 均 显 著 优 于 PSO 算法 和 ISM 算法 ， 


ISM 算法 和 PSO $ 
较 优 待 选 方案 。 
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